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Spin, charge, and orbital structures in models for doped manganites are studied by a combination 
of analytic mean- field and numerical relaxation techniques. At realistic values for the electron- 
phonon and antiferromagnetic ti g spin couplings, a competition between a ferromagnetic (FM) 
phase and a charge-orbital ordered (COO) insulating state is found for £=1/4, 3/8, and 1/2, as 
experimentally observed in Pri-^CaoiMnOs for x=0.3~0.5. The theoretical predictions for the spin- 
charge-orbital ordering pattern are compared with experiments. The FM-COO energy difference is 
surprisingly small for the densities studied, result compatible with the presence of a robust colossal- 
magnetoresistive effect in Pri-^Ca^MnOa in a large density interval. 
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The explanation of the colossal magnetoresistance 
(CMR) phenomenon in perovskite manganese oxides is 
currently one of the moat challenging problems in con- 
densed matter physics.u In addition to their unusual 
magneto-transport properties, manganites present a very 
rich phase diagram involving phases with spin, orbital, 
and charge order. Early investigations of these materials 
relied on the "double-exchange" (DE) mechanism for fer- 
romagnetism, in which e g electrons optimize their kinetic 
energy if the t 2g spins background is polarized. How- 
ever, this simple picture is not sufficient to rationalize the 
CMR effect since the properties of the insulating state in- 
volved in the metal-insulator transition play a key role. 
In fact, the understanding of the A-type antiferromag- 
netic (AFM) insulator LaMn03 requires a two orbital 
model and strong clcctron-phonon or Coulomb interac- 
tions to induce the complex spin and orbital arrangement 
characteristic of this state.H The evolution of the undoped 
phase with light hole doping, and its eventual transforma- 
tion into a ferromagnetic (FM) metal,0 is nontrivial and 
it is expected to proceed through a mixed phase process 
involving nanometer domainsD The conspicuous percola- 
tive characteristics of the transition have been recently 
explained by the influence of quenched disopder in the 
hopping and exchanges of manganite models B 

At high hole densities near x=l/2 a strong CMR phe- 
nomenon occursfl For example, in Ndi/2Sr 1 / 2 Mn03,Ll a 
low magnetic field of a few Teslas is enough to induce 
a metal insulator transition between the FM metallic 
and charge-ordered (CO) insulating phases. Based on re- 
cent computational studies at x—1/2, that clearly showed 
the presence of a first-order level-crossing transition be- 
tween the FM and CE-type CO statesja a simple picture 
to explain this result can be constructed. At densities 
slightly above the critical hole doping x c ~ 1/2 of the FM- 
CO transition it is reasonable to expect that small mag- 
netic fields will transform the ground state from CO to 
FM, since these states are close in energy. However, the 
FM metallic and CO insulating phases have a quadratic 



and linear dependence with x, respectively, and small 
magnetic fields will become rapidly ineffective to pro- 
duce such a transition as x grows away from x c . Thus, 
in this simple scenario the CMR effect can only occur 
in a narrow density window around x c . Unfortunately, 
some experiments suggest that this scenario is incom- 
plete. In_fact, it is well-known that for Pri-^Ca^MnOs 
(PCMO)El the CMR effect occurs in a wide density range 
0.3<x<0.5, and at first sight it appears that a simple 
level-crossing picture is not suitable for this compound. 

It is the purpose of this paper to discuss a possi- 
ble alternative explanation for the CMR phenomenon in 
PCMO. In the new scenario it is still claimed that at a 
fixed density small magnetic fields can induce FM-CO 
transitions as in the simple level crossing picture, but 
the key new idea is that different types of CO phases are 
stabilized at different densities. The calculations are here 
carried out at some special dopings, £=1/4, 3/8, and 1/2, 
for technical reasons to be discussed below, and using the 
two-orbital model strongly coupled to the Jahn- Teller 
(JT) distortions. A remarkable result observed in the 
present study is that, for realistic parameters, the region 
of competition between FM and CO phases is found to 
occur virtually independently of x, namely for the same 
values of couplings. The FM-CO energy difference is sur- 
prisingly small, as long as the CO state is optimized at 
each density. These novel candidate CO phases compet- 
ing with the FM state are here described in detail. In 
previous literature it is usual to find references to these 
states as "x=l/2 CO plus defects", due to their experi- 
mentally observed similarities with the x=l/2 charge ar- 
rangement, although they have different electronic den- 
sities. Here concrete examples of CO states with x^l/2 
are presented. The discussion below lead us to believe 
that similar conclusions would be obtained if Coulomb 
interactions were included in the calculations. 

Let us consider the hopping of e g electrons, tightly cou- 
pled to localized spins and the JT distortions of the 
MnOe octahedra. Their Hamiltonian is 
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—c\ ba Cibcr), where c\ a a {ciba) is the annihilation operator 
for an e g electron with spin er in the d x i_ y i (d 3z 2_ r 2) 
orbital at site i, a is the vector connecting nearest- 
neighbor sites, and , is the hopping amplitude be- 
tween 7- and 7'-orbitals in neighboring sites along the 
a-direction, given by t* a =~V3t* b =- V3t% a = 3tf b =3t/4 
for a=x, tZ a =V3t y ab =V3t y ba = 3^ b =3t/4 for a=y, and 
%=3t/4, tl=tl b =t* b =Q for a=z. The Hund coupling 
Jh(> 0) links the e g electron spin Si=J] 7 a/3 c i 7 a <T Q/3Ci 7 /3 
and the localized i 2g spin Si assumed classical (|Si|=l). 
J' is the AFM-coupling between nearest-neighbor i 2g 
spins. The dimensionless clcctron-phonon coupling con- 
stant is A. Q21 and Q^i are, respectively, the (x 2 —y 2 )- 
and (3z 2 — r 2 )-type JT modes of the Mn0 6 octahedron. 

To simplify the model the widely used limit Jh=oo is 
considered in this work. In such a limit, the e g electron 
spin perfectly aligns along the i 2g spin direction, reduc- 
ing the number of degrees of freedom. A is evaluated 
from X—^/2EjT/t, where r-the static JT energy Ejt is 
estimated as -EjT=0.25eVtLil and t is typically 0.2eV in 
narrow-band manganese oxides such as PCMO. Thus, in 
this paper A is fixed as 1.6 to compare the theoretical 
predictions with experiments, but the results shown be- 
low are not much sensitive to small changes in A.li3 In the 
rest of the paper, the importance of the ./'-dependence of 
ground state energies is emphasized to address the FM- 
CO competition. 

To solve the Hamiltonian (1), in this paper both a nu- 
merical relaxation technique and an analytic mean-field 
(MF) approximation are employed. In the former, the 
optimized JT distortion and i 2g spin configuration are de- 
termined numerically by the Simplex Method. Although 
the result is very accurate, considerable CPU time is 
needed to achieve convergence, and it is difficult to treat 
large clusters with this method. In fact, the cluster stud- 
ied here is a 4x4x4 cube with periodic boundary condi- 
tions (PBC). However, this lattice size is enough for the 
present investigation at £=1/4, 3/8, and 1/2. If other 
dopings such as x=0.3 and 0.4 are studied, larger size 
lattices should be treated due to the complexity of the 
CO states described here, which typically have a large 
unit cell. In the analytic approach, on the other hand, 
the JT distortion is determined self-consistently at each 
site using the relations Q 2 i=— A(r a; i) and Q3i=— A(r z i), 
where the bracket denotes the average value. In this ap- 
proach, energies for a variety of possible t 2g spin pat- 
terns are compared to find the lowest-energy state. This 
method has the advantage that a large-size lattice can be 
treated without much CPU time, but it must be checked 
whether the obtained structure indeed corresponds to the 
lowest-energy state. Such a check is here carried out by 



comparing MF and unbiased numerical results. Thus, the 
combination of the analytic and numerical techniques is 
powerful to obtain accurate predictions rapidly. In addi- 
tion, the reliability of the relaxation technique has been 
checked by comparing data with unbiased Monte Carlo 
(MC) simulations.^] At x=0, the agreement between the 
relaxation technique and MC simulation was excellent in 
any dimensions within errorbarsjj indicating the reliabil- 
ity of the present combined analytic-numeric approach. 
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FIG. 1. Energies per site vs. J' of several spin arrange- 
ments for (a) x=l/4, (b) 3/8, and (c) 1/2 in a 4x4x4 cubic 
lattice with PBC. The solid circles and lines are obtained 
by the relaxation method and the MF approximation, re- 
spectively. The meanings of the lines are as follows: Thick 
solid (FM), thick broken (pCE), thin broken (A-AFM), thick 
dashed (rCE), thin dashed (sCE), thick dot-dashed (CE), thin 
dot-dashed (C-AFM), and dotted (G-AFM). (d) Charge cor- 
relation function in the FM phase for 3;= 1/4 (solid circle), 3/8 
(open square), and 1/2 (solid diamond), (e) Schematic charge 
configuration for the FM insulating phase at 2=1/4. To save 
space, only a 2x2x4 cluster is shown, but the structure is 
periodically repeated in all directions, (f) Energy difference 
8E between the FM and CO phases vs. x for J'=0.09. The 
broken curve is a spline interpolation. Note that the region 
with positive SE corresponds to a CO ground state. 

Now let us analyze our results for x=l/4, 3/8, and 
1/2 (Fig.l(a)-(c)), studying energies vs. J' for several 
spin patterns. The MF results agree accurately with 
those of the relaxation technique, confirming the relia- 
bility of the MF method. For J'< 0.1, the FM phase 
is stabilized and it is insulating at a;=l/4, as shown 
in Fig. 1(d) where sharp peaks can be observed in the 
Fourier transform of the charge correlation function de- 
fined as n(q) = (l/A)X; i!j e" lcl -( i - j )((ni-n)(nj-n)). Here 

N is the total number of sites, rii=Y^ Ci 7<T , and 
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n(=l— x) is the average electron number per site. The 
peaks in Fig. 1(d) indicate a CO pattern, schematically 
shown in Fig. 1(e), in which two-dimensional (2D) planes 
witlya=l and 1/2, respectively, are stacked along the z- 
axis.113 This charge-ordered FM insulating phase emerg- 
ing from our calculations may be experimentally detected 
for PCMO. Note that the "FM insulating" phase of man- 
ganites has not been analyzed in detail in previous stud- 
ies, and here it is conjectured that it has charge-ordering. 
For £=3/8 and 1/2, the FM phases are found to be metal- 
lic, since no clear peak can be observed in n(q). This 
insulator-metal change in the FM state as a function of 
x agrees quite well with PCMO experiments. 

Around J'«0.1, a close competition occurs among the 
FM phases and several CE-type CO states (to be de- 
scribed later), indicating that a small perturbation can 
easily induce a first-order transition between FM and CO 
phases. It should be emphasized the remarkable result 
that such a competing region does not sensitively depend 
on x, suggesting that a CMR phenomenon can occur in a 
wide range of x as observed in PCMO for x = 0.3^0.5. In 
fact, as shown in Fig. 1(f), the energy difference between 
the FM and CO phases are 0.004 at x=3/8 and 0.018 at 
1=1/2 for J'=0.09. If t is assumed to be 0.2eV, those 
are about 8 and.-3j6 Tesla, similar values as observed in 
the experimentsJU 

Let us discuss the detail structure of the CO phase 
around J'~0.1. Our results suggest that the possible CO 
phases are given by several combinations of the 2D CE- 
type AFM spin pattern. To simplify the discussion, the 
CE-type AFM configuration in Fig. 2(a) is used as the 
basic pattern. The possible ground states are classified 
into four types by the stacking manner along the z-axis. 
(i)"the planar CE-type" (pCE), in which the pattern (a) 
stacks along the z-axis without any change, (ii) "CE- 
type" (CE), in which the patterns (a) and (b) stack along 
the z-axis alternatively. Note that this is the CE-type 
AFM structure observed in several half-doped mangan- 
ites. (iii) "shifted CE-type" (sCE), in which patterns 
(a) and (c) stack along the z-axis alternatively. Note 
that (c) is obtained by shifting (a) by one lattice spacing 
along the y-axis. This structure, with no charge-stacking, 
was suggested as a possible x—1/2 ground state, if the 
NaCl-type_CO occurs due to the strong nearest-neighbor 
repulsion.EJ (iv) "rotated CE-type" (rCE), in which pat- 
terns (a) and (d) stack along the z-axis alternatively. 
Note that (d) is obtained by rotating (a) by tt around 
a certain corner site in the zigzag FM chain in the x-y 
plane, or by a two lattice spacing shift. Both the rCE 
and sCE states have the same number of AF and FM 
links and their energies are J'-independent. 

As observed in Figs.l(a)-(c), both the pCE- and rCE- 
phases compete with the FM phase for 2=1/4 and 3/8, 
while at £=1/2, both sCE- and, especially, the CE-phase 
are competitive with the FM phase. As for the CO pat- 
tern, the CE-, pCE-, and rCE-phases exhibit the charge 
stacking due to a peak in n(q=(n 7 tt, 0)) (Fig. 2(e)). On 
the other hand, the sCE-phase has the NaCl-type charge 



ordering with a peak at q=(7r, tt, tt). Since experimen- 
tally z-axis charge stacking has been observed, the pCE-, 
rCE-, and CE-type phases are the best candidates for the 
ground state of the CO phases, since their energies are 
low in the analysis reported here. 




FIG. 2. (a)-(d) Four types of AFM CE-type spin configura- 
tions in the x-y plane. Solid and open symbols indicate up and 
down t2 g spins, respectively. The lobes indicates (3a; 2 — r 2 )- 
or (3j/ 2 — r 2 )-orbital at Mn 3+ site, while the circles denotes 
either Mn 4+ or an imperfectly doped site. The stacking along 
the z-axis of (a)- (a), (a)-(b), (a)-(c), and (a)-(d) lead to pCE-, 
CE-, sCE-, and rCE-phases, respectively, (e) Charge corre- 
lation for pCE-, CE-, sCE-, and rCE-phases. For pCE- and 
rCE-types, results at £=3/8 are shown, while for sCE- and 
CE-types, x is chosen as 1/2. (f) feg spin correlations for 
pCE-, CE-, sCE-, and rCE-phases. 

Now consider the spin correlations 5(q)=(l/iV) ^\ j 
e — *«!•(»— jj (Sj-Sj). The spin patterns of the states found 
by combining Fig. 2(a)-(d) have the characteristic zigzag 
FM chains in the x-y plane, leading to peaks at (tt, 0, q z ) 
and (tt/2, tt/2, q z ), where q z depends on the stacking ar- 
rangement. For pCE- and CE-types, q z is given by and 
tt, respectively, due to the FM- and AFM-spin stacking 
along the z-axis (Fig. 2(f)). For rCE-type, the peaks ap- 
pear at (tt, 0,0) and (tt/2, tt/2, tt), since half the z-axis 
bonds are FM and half AF. For sCE-type, four peaks 
exist at (77, 0, q z ) and (tt/2, 7r/2, q z ) with q z =0 and tt. hi 
neutron experiments for PCMO with x — 0.3 and 0.4,li3 
a peak was found at q=(7r/2, 7r/2, 0), suggesting that the 
pCE-state is the best candidate for the ground state in 
0.3 0.4. The well-defined charge and spin arrange- 

ment discussed here considerably improves over previous 
more vague "2;= 1/2 CE-type plus defects" descriptions of 
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this state. At x=l/2, the (n/2, 7r/2, tt) peakt3 indicates 
the CE-type as the ground state. 

Note that the present results do not always provide the 
pCE-type as the lowest-energy state in a sizable region of 
J' at £=1/4 and 3/8. To stabilize the pCE-phase it may 
be necessary to include the Coulomb interactions, espe- 
cially the nearest-neighbor repulsion V. In fact, if the 
nearest-neighbor charge correlation Cnn= X^<i j)(PiPj) is 
evaluated, Cnn for the pCE phase is found to be smaller 
than that of the rCE-type, indicating that rCE-type is 
indeed more energetically penalized by V. Thus, it is rea- 
sonable to expect that a stability "window" for the pCE- 
state will appear if it were possible to include accurately 
the Coulomb interactions. Note, however, that the size of 
such a "window" will sensitively depend on the parame- 
ter choice, cluster size, and approximations, since a small 
perturbation can easily modify ground state -properties in 
a region where several states are competing Jill Moreover, 
at £=1/2, the sCE-phase with the NaCl-type CO will be 
stabilized if V is very large, making the situation more 
complicated since it is known that charge stacking, not 
present in the sCE- phase, occurs at x=l/2. Nevertheless, 
from our present analysis it can be safely concluded that 
either the pCE- or rCE-states can be the lowest-energy 
state and play the role of the experimentally observed 
CO phase in the interesting region of PCMO. 

As for the orbital ordering, the alternation of 
(3x 2 — r 2 )- and (3y 2 — r 2 )-orbitals appears in the x-y 
plane, although the exact orbital shape is slightly de- 
pendent on the stacking manner. At the corner sites 
in the FM zigzag path, except for sCE-phase, our cal- 
culation suggests that the orbital is polarized along the 
z-axis and the (3z 2 — r 2 )-orbital is occupied—This result 
agrees well with the experimental data.EjO To distin- 
guish between pCE- and rCE-phases, the stacking of or- 
bitals along the z-axis is important. Namely, for the pCE 
and rCE phases, the (3x 2 —r 2 )- and (3y 2 — r 2 )-orbitals 
stacks in ferro and antiferro patterns, respectively. This 
difference can be detected for PCMO at £~3/8 by using 
the resonant X-ray scattering technique. 

Finally, the possibility of spin canting is here briefly 
discussed. For pCE-type at £=3/8, the t2 S spins cant- 
ing between two adjacent x-y planes along the z-axis was 
studied as a possible way to convert continuously from 
pCE-type to CE-type. However, when the ground state 
energy is plotted against the canting angle, local minima 
were found only at angles corresponding to the extreme 
cases pCE- and CE-types, and the spin canting phase 
was not stabilized in the present work. However, the 
existence of the spin canting state cannot be totally ex- 
cluded if a finite value of Jh is considered (here Jh=oo 
was used). 

Summarizing, spin-charge-orbital ordering in mangan- 
ites has been investigated at £=1/4, 3/8, and 1/2 using 
analytical and numerical techniques. For fixed values of 
A and J', it has been found that the energy difference 



between the FM and CO phases is remarkably small for 
£=1/4, 3/8, and 1/2, and a small magnetic field can in- 
duce a CO-FM transition. This result is a first step to- 
ward a possible explanation of the CMR effect in PCMO 
which occurs in a robust density range. Contrary to pre- 
vious descriptions of these states as "£=1/2 plus defects" , 
here it has been shown that specific charge arrangements 
without defects can be constructed to represent the CO 
states at x^O. Our predictions can be verified using X 
ray scattering experiments. 
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